home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / CGBSL.z / CGBSL
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. CCCCGGGGBBBBSSSSLLLL((((3333FFFF))))                                                            CCCCGGGGBBBBSSSSLLLL((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CGBSL   - CGBSL solves the complex band system A * X = B  or  CTRANS(A) *
  10.      X = B using the factors computed by CGBCO or CGBFA.
  11.  
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  14.       SUBROUTINE CGBSL(ABD,LDA,N,ML,MU,IPVT,B,JOB)
  15.  
  16. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  17.      On Entry
  18.  
  19.      AAAABBBBDDDD COMPLEX(LDA, N)
  20.         the output from CGBCo or CGBFA.
  21.  
  22.      LLLLDDDDAAAA INTEGER
  23.         the leading dimension of the array  ABD .
  24.  
  25.      NNNN INTEGER
  26.         the order of the original matrix.
  27.  
  28.      MMMMLLLL INTEGER
  29.         number of diagonals below the main diagonal.
  30.  
  31.      MMMMUUUU INTEGER
  32.         number of diagonals above the main diagonal.
  33.  
  34.      IIIIPPPPVVVVTTTT INTEGER(N)
  35.         the pivot vector from CGBCO or CGBFA.
  36.  
  37.      BBBB COMPLEX(N)
  38.         the right hand side vector.
  39.  
  40.      JJJJOOOOBBBB INTEGER
  41.         = 0         to solve  A*X = B ,
  42.         = nonzero   to solve  CTRANS(A)*X = B , where
  43.         CTRANS(A)  is the conjugate transpose.  On Return
  44.  
  45.      BBBB the solution vector  X . Error Condition
  46.  
  47.      AAAA division by zero will occur if the input factor contains a zero on the
  48.      diagonal.  Technically this indicates singularity but it is often caused
  49.      by improper arguments or improper setting of LDA .  It will not occur if
  50.      the subroutines are called correctly and if CGBCO has set RCOND .GT. 0.0
  51.      or CGBFA has set INFO .EQ. 0 .  To compute  INVERSE(A) * C  where  C  is
  52.      a matrix with  P  columns
  53.         CALL CGBCO(ABD,LDA,N,ML,MU,IPVT,RCOND,Z)
  54.         IF (RCOND is too small) GO TO ...
  55.         DO 10 J = 1, P
  56.         CALL CGBSL(ABD,LDA,N,ML,MU,IPVT,C(1,J),0) 10 CONTINUE LINPACK.  This
  57.      version dated 08/14/78 .  Cleve Moler, University of New Mexico, Argonne
  58.      National Lab.  Subroutines and Functions BLAS CAXPY,CDOTC Fortran
  59.      CONJG,MIN0
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.